Parity Nonconservation in the Photodisintegration of the Deuteron at Low Energy 
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The parity-nonconserving asymmetry in the deuteron photodisintegration, 7 + d — > n + p, is 
considered with the photon energy ranged up to 10 MeV above the threshold. The aim is to 
improve upon a schematic estimate assuming the absence of tensor as well as spin-orbit forces in the 
nucleon-nucleon interaction. The major contributions are due to the vector-meson exchanges, and 
the strong suppression of the pion-exchange contribution is confirmed. A simple argument, going 
beyond the observation of an algebraic cancellation, is presented. Contributions of meson-exchange 
currents are also considered, but found to be less significant. 



I. INTRODUCTION 

Some interest, both experimental and theoretical, has 
recently been shown for the study of parity nonconser- 
vation in the deuteron photodisintegration by polarized 
light. Historically, it was its inverse counterpart: the 
net polarization in radiative thermal neutron capture by 
proton, n + p — > d + 7, which attracted the first at- 
tention The experimental study was performed by 
the Leningrad group, taking advantage of new techniques 
measuring an integrated current 0. The non-zero polar- 



ization obtained, P~, 



-(1.3 ±0.45) x lCT 6 , motivated 



many theoretical calculations in the frame of strong and 
weak interaction models known in the 70's (see for in- 
stance Refs. 0,0, 01)- The theoretical results were con- 
sistently within the range P 1 = (2 ~ 5) x 10~ 8 , which is 
smaller than the measurement by a factor of 30 or more 
in magnitude and, moreover, of opposite sign. The diffi- 
culty to understand the measurement and, also perhaps, 
the novelty of the techniques, which have been exten- 
sively used later on, led to a special reference to this 
work as "Lobashov experiment". 

Later estimates with modern nucleon-nucleon (NN) 
potentials, both parity-conserving (PC) and parity- 
nonconserving (PNC), give values of P 1 roughly within 
the same theoretical range as above. On the experimen- 
tal side, new results were reported in the early 80's by 
the same Leningrad group ^giving P 7 < 5 x 10~ 



P 7 = (1.8 ± 1.8) x 10" 



j| and 

Practically, these results 
indicate an upper limit of P 7 , which is not very constric- 
tive. Since Leningrad group's last report, the "Lobashov 
experiment" has long been forgotten by both experimen- 
talists and theorists. Recent experiments such as elas- 
tic p-p scattering (TRIUMF 0) and polarized thermal 
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neutron capture by proton (LANSCE 0), which directly 
address the problem of PNC NN interactions; and quasi- 
elastic e-d scattering (MIT-Bates which indi- 
rectly involves these interactions, have however raised a 
new interest for the study of PNC effects in few-body 
systems. In what could be a golden age for these studies, 
the "Lobashov experiment" is again evoked. 

While it seems that there is not much prospect for per- 
forming the "Lobashov experiment" in a near future, the 
inverse process, on the contrary, could be more promis- 
ing. In this reaction, 7 + d — > n + p, where a deuteron is 
disintegrated by absorbing a circularly polarized photon, 
it is expected that, near threshold, the PNC asymmetry 
(Ay) is equal to the polarization in the "Lobashov exper- 
iment". This last one can thus be tested from a different 
approach. 

The asymmetry Ay in the deuteron photodisintegra- 
tion was first calculated by Lee ^3] up to the photon 
energy w 7 ~ 3.22 MeV, which is 1 MeV above the thresh- 
old. In this energy domain, where the dominant regular 
transition is Ml, the result was within the theoretical 
range of P 7 . Later on, Oka extended Lee's work, up to 
ui 1 ~ 35 MeV Though the cross section still re- 

ceives a contribution from the Ml transition, the dom- 
inant contribution comes from the El transitions. This 
offers a pattern of PNC effects different from the one at 
very low photon energy. It was found that Ay shows a 
great enhancement at w 7 > 5 MeV, mainly due to the 
PNC 7r-exchange contribution. If such an enhancement 
were observed in the experiment, it would provide an 
important and unambiguous determination of the weak 
irNN coupling constant h\. However, a recent schematic 
calculation of Ay by Khriplovich and Korkin 0], partly 
suggested by one of the present author, showed critical 
contradiction to Oka's result, with a huge suppression of 
Ay at the energies u 1 > 3 MeV. 

On the experimental side, a measurement of the asym- 
metry Ay in 7 + d — > n+p was considered in the 80 's by 
E. D. Earle et al. 0,E3] but no sensitive result was re- 
ported. However, due to advances in experimental tech- 
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niques and instrumentation, the measurement of A n be- 
comes more feasible nowadays and several groups at JLab 
H3, IAS A (Athens), LEGS (BNL), TUNL, and SPring-8 
show interest in such a measurement. It is therefore im- 
portant to understand and improve previous estimates. 

In this work, we carefully re-examine the 7 + d — > n+p 
process with two main purposes: 

1. Determine how the enhancement of the h\ contri- 
bution in Oka's results will change when the cal- 
culation is completed with missing parity-admixed 
components in the final state, in particular in the 
3 Pi channel. The role of this last one was revealed 
by the schematic estimate of Ref. 0] . 

2. Determine the uncertainty of Khriplovich and Ko- 
rkin's calculation in which very simple wave func- 
tions are used. 

It is straightforward to deal with the point 1. In Ref. 
[lH |. a nice and simple argument about the cancella- 
tion of the h\ contribution from the final 3 Po, 3 Pi, and 
3 p2 states along with their parity-admixed partners was 
given. However, the argument assumed the absence of 
tensor as well as spin-orbit forces, which are important 
components of the AW interaction. In order to address 
these two points (missing components and simplicity of 
the wave functions), we elaborate our calculation with 
the Argonne Vis NN interaction model. We thus include 
the 1 Sq, 3 Po, 3 Pi, 3 P2- 3 F 2 channels, deuteron D-state, 
and all their parity-admixed partners consistently. They 
represent a minimal set of states that allows one to verify 
the results of the schematic model as well as to include 
the effect of the tensor and spin-orbit forces that mani- 
fest differently in these various channels. We also include 
other channels, whose role is less important however. As 
for the El operator, we employ the Siegert's theorem 
[lfi| . which takes into account the contribution of some 
PC and PNC two-body currents. The small photon en- 
ergy considered here (w 7 < 12 MeV) justifies this usage. 
Since there is no theorem similar to the Siegert one for 
the Ml transition operator, two-body currents have to 
be considered explicitly for both the PC and PNC parts. 
Adopting Desplanques, Donoghue, Holstein (DDH) po- 
tential of the weak interaction the asymmetry A~ t 
will be expressed in terms of the weak irNN, pNN and 
ujNN coupling constants, with corresponding coefficients 
indicating their relative importance. 

This paper is organized as follows. In Sect. [H] we 
review the basic formalism underlying the calculation, 
which involves both one- and two-body currents. In Sect. 
IIHI we show the results and some discussions follow. A 
particular attention is given to a comparison with earlier 
works and to new contributions from PNC two-body cur- 
rents. A simple argument explaining the suppression of 
the pion-exchange contribution is also given. Conclusions 
are given in Sect. IIVI An appendix contains expressions 



of El and Ml transition amplitudes due to the PNC 
two-body currents considered in the present work. 
II. FORMALISM 

For a photodisintegration of an unpolarized target, the 
asymmetry factor is defined as 



-4, 



where &+(-) denotes the total cross section using right- 
(left-) handed polarized light. By spherical multipole ex- 
pansion, it could be expressed as 



i4.™ 



2 Re ^ \f ej F M j 5 + Fmj F E j 5 
f,j,J 

f,i,J 



(1) 



In this formula, the normal electromagnetic (EM) and 
PNC-induced EM form factors, F X j and F X j 5 , with X 
and J denoting the type and multipolarity of the transi- 
tion between a specific initial (i) and final (/) states, are 
defined in the same way as Refs. USE!. They depend 
on the momentum transfer q, which equals to the photon 
energy uo 1 in this current case. The form factors F X j 5 
(and so does the asymmetry) vanish unless some PNC 
mechanism induces parity admixtures of wave functions 
and axial- vector currents. 

In this work, we consider the photon energy lu 1 = q 
up to 10 MeV above the threshold. As the long wave- 
length limit, (qr) C 1, is a good approximation, the 
inclusion of only dipole transitions, i.e. El and Ml, is 
sufficient. This leads to 10 possible exit channels con- 
nected to the deuteron state by angular momentum con- 
siderations. Among them, 1 So, via the Ml transition, 
and 3 Po, 3 Pi, 3 P2- 3 F2, via the El transitions, dominate 
the cross section. 

The transverse multipole operators assume a full 
knowledge of nuclear currents. This requires, besides the 
one-body current from individual nucleons, a com- 
plete set of two-body exchange currents (ECs) which 
is consistent with the nucleon-nucleon (NN) potential. 
These ECs are usually the sources of theoretical uncer- 
tainties, because the NN dynamics is still not fully un- 
derstood. While there is no alternative for the evaluation 
of Fmj, the Siegert theorem does allow one to trans- 
form the evaluation of Fej into the one of charge mul- 
tipole Fcj. The fact that the PC NN interaction does 
not give rise to exchange charges at 0(1) removes most of 
the uncertainties related to exchange effects: knowledge 
of the one-body charge is sufficient for a calculation 
good to the order of l/m N . 

In the framework of impulse approximation and using 
the Siegert theorem, one gets, for the deuteron photodis- 
integration (Ef — Ei — lo 7 = q and J, = 1), 
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where e, = e (1 + t/)/2 and /tj = e (/a s + n v rf)/2 with 
/x s = 0.88 and /U v = 4.70; Y and Y" are the spheri- 
cal and vector spherical harmonics. In these expressions, 
the approximated results are obtained by replacing the 
spherical Bessel function ji (q x) with its asymptotic form 
as q — > 0, i.e. qx/3, at the long wavelength limit and 
keeping terms linear in q (the lowest order); they could 
be related to the forms of (£1 (1 )) and (Ml^) often 
adopted in the literature. In our numerical calculation, 
the identity relations are employed instead. Note that 
the one-body spin current is conserved by itself and not 



constrained by current conservation. In Eq. (2J, this 
one-body spin current (2nd line) is of higher order in q 
compared with the Siegert term (1st line), however, it 
is kept for completeness. As for the PNC-induced form 
factors and , one only has to replace either 

the initial or final state by its opposite-parity admixture, 
(J f | or | Ji), and add a factor u i" for El or 11 -i" for Ml 
matrix elements (in relation with our conventions). 

The non-vanishing matrix elements for the five domi- 
nant exit channels are thus 



J 



1. ^q: 



(MlW) 
(El?) 



fJ-v 



drU*( 1 S )U d ( 3 S 1 ) 



rdrU*( 3 P ) \u d ( 3 S 1 )~V2U d ( 3 D 1 



V3 



rdrU*CS )U d ( 1 P 1 ) 



(4) 



(5) 



2. 3 P : 



(#) 
(Ml«) 



rdrU*( 3 P ) f£/ d ( 3 Si)-\/2[/ d ( 3 £>i) 



dr 



U*CS ) ?7 rf ( 3 5i) - -j= U*( 3 P ) U d ?Pi) 



j2f*s_J/2 i drU * { 3 Po) u d( 3 Pl) 

O 777, at 



(C) 



(7) 



3. 3 Pi: 



(£1«) 



(Ml™) 



7=3 J ' rdrU * {3Pl) 



U d ( 3 Si) + -j=U d (! i D 1 ) 



-i^L f drU *{ 3 P 1 )U d Q;P{)-i^±±ll / drU*( 3 P 1 )U d ( 3 P 1 ) 
m N J V2m N 



— i ■ 



m 



rfrt/*( 3 5i)l/ d ( 3 ^i)+ ^ s / _ 3/2 / drU*( 3 D 1 )U d ( 3 D 1 ). 
V2m JV 



(8) 
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4. 3 P 2 - 3 F 2 



(El 



(Ml 



W\ = 



V5 



dr\ U*{ 6 P 2 ) 



U d { 3 Sy)-jj=U d { 3 D{) 



+^U*( 3 F 2 )U d ( 3 D 1 ) 



(10) 




U*( 3 P 2 ) UdCPi) - ^U*CD 2 ) U d ( 3 D 1 ) 



U*( 3 P 2 ) U d ( 3 P 1 ) + [Tf £> 2 ) [7,(3^) 



(11) 



r 



Results for the remaining five less important channels 
( 3 S , i- 3 Di, 1 Pi, 1 D 2 , 3 £>2) will be included in numerical 
works. The r- weighted radial wave functions for scatter- 
ing and deuteron states, U and U d) along with their par- 
ity admixtures, U and U d , are obtained by solving the 
Schrodinger equations. Details could be found in Ref. 

m- 

By taking the square of normal EM form factors (PC 
response function) or the product of normal and PNC- 
induced ones (PNC response function), we can directly 
compare Eqs. (5a-5h) in Ref. 0|. After removing fac- 
tors due to wave-function normalizations, the differences 
are: 

1. The parity admixture of the scattering 3 Pi state is 
included in our work: The admixtures U( 3 Si) and 
U( 3 Di) are solved from the inhomogeneous differ- 
ential equations with the source term modulated by 
U( 3 Pi). They are not orthogonal to the deuteron 
state and thus should not be ignored. Actually, 
they are required to ensure the orthogonality of the 
deuteron and the 3 P\ scattering states once these 
ones are allowed to contain a parity-nonconserving 
component. 

2. The terms involving the scalar magnetic moment 
are different: Looking for instance at the Ml ma- 
trix element between U( 3 Pi) and U d ( 3 Pi), the ef- 
fective Ml operator is proportional to /x s S + L/2. 
By the projection theorem, (S) = (L), the overall 
factor should be /i s + 1/2, not /i s + 1 as in Ref. 

R looks as if this work ignored the 1/2 fac- 
tor in front of the L operator. 

Both points involve the spin-conserving PNC interaction, 
which is dominated by the pion exchange. Therefore, how 
these differences change the sensitivity of A 1 with respect 
to h\ will be elaborated in next section. 

Now we discuss, in two steps, extra contributions due 
to ECs when one tries to go beyond the impulse approx- 
imation together with the Siegert-theorem framework. 

First, when PC ECs are included, their contribution 
to Ml matrix elements, F^[, definitely needs to be cal- 
culated. On the other hand, as PC exchange charges are 



(51 

higher-order in the nonrelativistic limit, F^ t is supposed 
to take care of most two-body effects, and the remaining 

(2) 

contribution AFg/ can be safely ignored. This argument 
also applies for the PNC-induced form factors involving 

r (2) 

the PC ECs: one needs to consider F^i but can leave 
out Aif> 5 . 

Second, the inclusion of PNC ECs, to the first order 
in weak interaction, only affects the PNC-induced form 

factors. The contribution F^i is calculated by using 
the Ml operator constructed from the PNC ECs and 
unperturbed wave functions (so we use a prime to remind 
the difference from parity-admixture contributions) . One 
special feature of PNC ECs is that they do have exchange 
charges of O(l) H^]. Therefore, one should include them 

As a last remark, we note one advantage of nuclear 
PNC experiments in processes like photodisintegration 
or radiative capture. The real photon is "blind" to the 
nucleon anapole moment, which could contribute other- 
wise to PNC observables in virtual photon processes. Be- 
cause this P-odd T-even nucleon moment is still poorly 
constrained both theoretically and experimentally, the 
interpretation of real-photon processes, like the one con- 
sidered here, is thus comparatively easier. 



III. RESULTS AND DISCUSSIONS 

For practical purposes, we use the Argonne Vis jl^ 
(Auig) and DDH potentials as the PC and PNC 
AW interactions, respectively. In comparison with earlier 
works in the 70's or the 80's, a strong interaction model 
like Avis offers the advantage that the singlet-scattering 
length is correctly reproduced, due to its charge depen- 
dence. Correcting results with this respect is therefore 
unnecessary. 

The total cross section is plotted in Fig. as a func- 
tion of the photon energy and labeled as "IA+Sieg". Its 
separate contributions from El and Ml transitions are 
also shown on the same plot (labeled accordingly). The 
Ml transition only dominates near the threshold region; 
as the photon energy reaches about 1 MeV above the 
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Figure 1: The total cross section as a function of the photon 
energy. The main result is the curve labeled as "IA+Sieg", 
and the curves "El" and "Ml" showing contributions from 
corresponding transitions. The curve "IA" is the result of a 
pure impulse approximation calculation, where no two-body 
contribution is included. 



threshold, the El transition overwhelms. Away from the 
threshold, the calculated results agree well with both ex- 
periment and existing potential-model calculations up to 
10 MeV Hi)]- Such a good agreement shows the useful- 
ness of the Siegert theorem, by which most of the two- 
body effects are included. Compared with the curve la- 
beled by "IA", the result of impulse approximation, one 
sees the increasing importance of these two-body con- 
tributions as w 7 gets larger. On the contrary, because 
Ml matrix elements are purely one-body, we expect our 
near-threshold results smaller than experiment by about 
10% US]- This discrepancy, originally found in the radia- 
tive capture of thermal neutron by proton (the inverse of 
deuteron photodisintegration) , requires various physics 
such as exchange currents and isobar configurations, to 
be fully explained. Here, we qualitatively estimate a 5% 
error for the calculation of Fmi near threshold. 

When calculating the PNC-induced matrix elements 
with the DDH potential, we use the strong meson-nucleon 
coupling constants: g-nNN = 13.45, g pNN = 2.79, and 
9ujnn = 8.37, and meson masses (in units of MeV): = 
139.57, m p = 770.00, and = 781.94. The resulting 
asymmetry is then expressed in terms of six PNC meson- 
nucleon coupling constants h's as 

A, = ci hi + c 2 h° p + c 3 h\ + c 4 h 2 p + c 5 hi +c 6 hi, (12) 

where the six energy-dependent coefficients ci...g show 
the sensitivity to each corresponding coupling. It 
turns out that, for the energy range considered here, 
C2, C4, C5 ^> ci ^ C3, cq. This implies the asymmetry 
has a larger sensitivity to the isoscalar and isotensor cou- 
plings than to the isovector ones. The detailed energy 
dependences of these "large" and "small" coefficients are 
shown in Fig. [2 
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Figure 2: The energy-dependences of "large" coefficients C2, 
c 4 , and c 5 (top panel) and "small" coefficients ci, 03 , and ce 
(bottom panel) in the asymmetry parametrization, Eq. 112t . 



In principle, these results are independent. In prac- 
tice however, they can be shown to depend on three 
quantities, reflecting the dominant role of the various 
S *-> P neutron-proton transition amplitudes at low en- 
ergy. These amplitudes have some energy dependence 
which is essentially determined by the best known long- 
range properties of strong interaction models. They can 
therefore be parametrized by their values at zero energy 
IH H3- HH; including at the deuteron pole. To a large 
extent, they can be used independently of the underlying 
strong interaction model, quite in the spirit of effective- 
field theories that they anticipated [13] ■ In the case of 
the Avis model employed here, they are given by: 



m N A t = 
m N A s = 

m N C = 1.023 hi 



0.043 h° p - 0.022/i° 



0.125 h° - 0.109 hi + 0.102 h 2 p , 

0.007 h p - 0.021 hi. (13) 



The largest corrections to the above approach occur for 
the PNC pion-exchange interaction which, due to its long 



range, produces some extra energy dependence and siz- 
able P «-> D transition amplitudes. They can show up 
when the contribution of the S *-> P transition amplitude 
is suppressed, like in this work. 

For w 7 = 2.235 MeV, which is very close to the disin- 
tegration threshold, we get the asymmetry 

[-8.44 h° - 17.65 h 2 p + 3.63 

+0(c 1 ,c 3 ,c 6 )} x 10" 3 . (14) 

Using the DDH "best" values as an estimate, we got 
2.53 x 10 -8 . By detailed balancing, one expects 



A (t h) 



Ath) 



p(th) 

7 



that A$f h ' equals the circular polarization P~f h ^ observed 
in the radiative thermal neutron capture by proton, given 
the same kinematics. Though our result does not ex- 
actly correspond to the same kinematics as the inverse 
process usually considered (the kinetic energy of thermal 
neutrons ~0.025 eV), it agrees both in sign and order 

of magnitude with existing calculations of 

We also performed a similar calculation for the latter case 

with Avis, and the result is 

[-8.75 h° - 17.47 h 2 p + 3.39 h° 

+0{c 1 ,c 3 ,c 6 )] x HT 3 . (15) 

This is very close to the result of A 1 quoted above. 
It is noticed that our expression of A^ h) at very low 

energy, and therefore that one for p!f h \ contains a con- 
tribution from the one-pion exchange (see the low-energy 
part of the c\ coefficient given in Fig. EJ. This feature, 
which apparently contradicts the statement often made 
in the past that this contribution is absent in Py , is 
due to the incorporation in our work of the spin term in 
Eq. (0, which represents a higher order term in q. This 
correction also explains the difference in the behavior of 
the c\ coefficient with the Oka's result (T^ |. 

We note that, because the Ml transition dominates at 
the threshold and we only use the impulse approximation 
for its matrix element, there should be approximately 
a —5% correction to (also .Ry ) when two-body 

effects are included in Fmi ■ On the other hand, as Fei b is 
calculated using the Siegert theorem, it should be reliable 

up to the correction of Fg 1 from the PNC exchange 
charge at 0(1). 

When the photon energy gets larger, one can see im- 
mediately that the asymmetry gets smaller. A prediction 
using the DDH best values is shown in Fig. In this 
figure, as soon as the photon energy reaches 1 MeV above 
the threshold, the asymmetry drops by one order of mag- 
nitude. Moreover, the sign changes around oj 1 = 4 MeV. 
This implies that a higher sensitivity (~ 10 -9 ) is needed 
for any experiment targeting at the kinematic range away 
from the threshold. Our calculation is consistent with the 
work by Khriplovich and Korkin ^4[; but is widely dif- 
ferent from the one by Oka 0. In the following, we 
make a closer comparison with these works and, then, 
present results for the contribution of various PNC two- 
body currents considered for the first time. 
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Figure 3: The asymmetry by using the DDH best values. 



A. Comparison with Oka's work 

The major difference comes from the pion sector. In Ref. 
[T3 |. where the scattering wave functions are obtained 
from the zero-range approximation and the deuteron is 
purely a 3 Si state, a simple angular momentum consider- 
ation leads to a null contribution from pions. Our result 
shows that the more complex nuclear dynamics has only 
small corrections, so the asymmetry is not sensitive to 
h\. However, it is not the case at all in Ref. 0|: the 
pion exchange dominates the asymmetry with the coeffi- 
cient ci being one or two orders of magnitude larger than 
our result. 

This discrepancy could be illustrated by considering a 
case where w 7 is 10 MeV above the threshold. In the 
central column of Table Q] we list the pertinent PNC re- 
sponses due to the pion exchange among the 5 dominant 
exit channels. In the right column, we simulate what the 
outcome will be if the analytical results of Eqs. (5a-5g) 
in Ref. 0] are used, i.e. with different factors involving 
jjL s and no parity admixture of 3 Pi state as mentioned in 
Sec. |nj Comparing the totals from both columns, one 
immediately observes the simulated result is bigger by an 
order of magnitude. More inspection shows that, while 
the changes of the /j, s factors do alter each response some- 
what, the major difference depends on whether the big 
cancellation from the 3 Pi admixture is included or not. 
By adding contributions from other sub-leading channels, 
the total will be further downed by a factor of 2.5. Thus 
the overall difference is about a factor of 30. 



Comparison with Khriplovich and Korkin's 
work 



The vanishing of the 7r-exchange contribution in 
Khriplovich and Korkin's work [Tj| supposes that the 
El transitions from the deuteron state to the different 
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Transitions 


Eqs. PHTTll 


Eqs. (5a-5h) in Ref. [131 


3 Po^f> 


0.449 


-0.142 


3 Pi <-> P 


-3.217 


-4.383 


3 Pi © 


3.942 


not considered 


3 P 2 ^P 


-1.231 


0.389 


3 P2 <-> P 


-0.142 


0.045 


s F 2 ^t> 


-0.151 


0.048 


3 P 2 <-> P 


-0.019 


0.006 


Total 


-0.371 


-4.037 



Table I: The dominant PNC responses due to the pion ex- 
change for lj 7 10 MeV above the threshold (in units of 
10~ J x h^). The central column is calculated by Eqs . 
111} . while the right column by Eqs. (5a-5h) in Ref. [lj- The 
symbol P denotes the deuteron state. 



scattering states, 3 Po, 3 P\ and 3 P2, are the same, which 
implies that one neglects both the tensor and spin-orbit 
components of the strong interaction. As these parts of 
the force have large effects in some cases, it is important 
to determine how the above vanishing is affected when a 
more realistic description of the interaction is used. 

We first notice that the isoscalar magnetic operator, 
/x s S + L/2, can be written as /i s J + (1/2 — /i s )L. As 
the operator J conserves the total angular momentum, 
it follows that the El transitions from the deuteron state 
to the 3 Pq and 3 P2 states will be proportional to fx s — 
1/2, in agreement with Eqs. Q and (11 H . A similar 
result holds for the 3 Pi state. For this transition, one 
has to take into account that the J operator connects 
states that are orthogonal to each other, including the 
case where they contain some parity admixture. This 
unusual but interesting result was originally suggested 
by a similar result obtained by Khriplovich and Korkin 
for the 1 Sq and 3 Po states [lj|. They used it later on 
for the 7r-exchange contribution on the suggestion of one 
the present authors. Taking this property into account, 
one can check that the different ^-dependent terms in 
Eq. © combine so that the quantity, ^ s — 1/2, can 
be factored out. This explains the cancellation of the 
two largest contributions in Table Q] 3.942 and —3.217, 
approximately proportional to 2/x s = 1.76 and — (fx s + 
1/2) = -1.38. 

Further cancellation is obtained when one considers the 
sum of the 7r-exchange contributions to the asymmetry 
A 1 corresponding to the different P states. Taking into 
account the remark made in the previous paragraph, it 
can be checked that contributions from Eqs. {JJ, and 
(T1H are proportional to 2, 3 and —5 and 4, —3, and —1 
for the 3 Si and 3 D\ deuteron components respectively 
(assuming that the 3 P wave functions are the same). As 
can be seen in Tabled the dominant contributions, 0.449, 
0.725 (= 3.942 - 3.217) and -1.231 are not far from the 
relative ratios 2, 3 and —5, expected for the 3 Si deuteron 
component. Possible departures can be ascribed in first 
place to the 3 Di deuteron component. 



The above cancellation calls for an explanation deeper 
than the one consisting in the verification that the al- 
gebraic sum of different contributions cancels. An ar- 
gument could be the following. In the conditions where 
the cancellation takes place (same interaction in the 3 P 
states in particular), a closure approximation involving 
spin and angular orbital momentum degrees of freedom 
can be used to simplify the writing of the PNC part of 
the response function that appears at the numerator of 
Eq. (QJ. Keeping only the factors of interest here, the 
interference term of El and Ml matrix elements can be 
successively transformed as follows 



SR oc r l (Ms - |) V - <f <f ) | J, 

M 

oc Y, [( 3 Si\U d ( 3 S 1 ) + 



Ud( 3 Di) 
V2 



3(S-fY-S z 



x r l U (S ij -ctq j ) S-f\ 3 Si) 



oc Tr 
= 0. 



E/d( 3 Si) 



Ud( 3 Di 
V2 



(3(S-f) 2 -S 2 



S-f 
(16) 



The first line stems from retaining the isoscalar part of 
the magnetic operator proportional to (/i g — 1/2)L (it is 
reminded that the J part does not contribute) . The next 
line is obtained by expressing the PC and PNC parts of 
the deuteron wave function as some operator acting on 
a pure \ 3 S\) state. Once this transformation is made, it 
is possible to replace the summation over the deuteron 
angular momentum components, M, by the spin ones, 
m s , which is accounted for at the third line. The last 
line then follows from the fact that the trace of the spin 
operator, S, possibly combined with a AS = 2 one, van- 
ishes. A result similar to the above one can be obtained 
for some contributions involving MECs. It is however no- 
ticed that some corrections involving the spin-orbit force, 
or spin-dependent terms in the El transition operator, 
which both contain an extra S factor in the above equa- 
tion, could lead to a non-zero trace and therefore to a 
relatively large correction. Of course, the above cancel- 
lation relies on the fact that no polarization of the initial 
or final state is considered. Had we looked at an observ- 
able involving such a polarization, like the asymmetry 
in the capture of polarized thermal neutrons by protons, 
the result will be quite different. As is well known, this 
observable is dominated by the 7r-exchange contribution 
01- 



C. Contributions of PNC ECs 

In Section El the contributions of PNC ECs were sum- 
marized in two additional PNC-induced form factors, 



b(S') 



F 



-Els 



and F 



(2') 
Ml 5 ' 



Now, we estimate these contributions 



by considering only the dominant channels 1 S'o, 3 Po, 3 Pi 
and 3 P2- 3 P2- As PI5 connects states of same parity, 
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only 1 Sq is allowed; therefore, plays a more impor- 

tant role for Ay near the threshold. On the other hand, 
MI5 connects states of opposite parity, which requires 
the other four channels, so F y M ^ has more impact on A 1 
at higher energies. The full set of PNC ECs which is con- 
sistent with the DDH potential was derived in j^, Eqs. 
(17-24). The whole evaluation is straightforward, how- 
ever tedious, so we defer all the analytical expressions in 
Appendix B and only quote the numerical results here. 

With the same parametrization as Eq. 111211 . the addi- 
tional contributions to the asymmetry by PNC ECs, via 
EI5 and MI5 respectively, are 

A 7 (<>) = ef^ + cf^, 

+c^hl + cphl. (18) 

The detailed energy-dependence of each coefficient is 
shown in Fig. 0] 

~(S') ~(2') 

The dominance of F^u near t ^ le threshold and F M ^ 
at higher energies could be readily observed in these 
plots. We discuss their significances to the total asym- 
metry in the following. 

For the case where the photon energy is 0.01 MeV 
above the threshold, only c 2 and c\ are substantial. 
The former coefficient is about 20% of C2, while the lat- 
ter one is only 2% of C4. By using the DDH best values, 
these contributions give an asymmetry about 1.4xl0~ 9 , 
which is a 6% correction. This is typically the order of 
magnitude one could expect from the exchange effects. 

is') 

As the energy gets larger, while the coefficients c 2 

( s') 

and c\ keep stable, the coefficients associated with 
A/I5 matrix elements grow linearly, roughly. The 
fastest growing one is c\ because the long-ranged pion- 

exchange dominates the matrix elements. Compara- 

(2') 

tively, c 2 has a smaller slope due to less overlap be- 
tween the effective ranges pertinent to the deuteron wave 
function and the p-exchange. 

For the case where the photon energy is 10 MeV above 

(2') (2 f ) (S') 

the threshold, c\ , c 2 , and c 2 are substantial. The 
first coefficient is about 50% of ci , and the latter two com- 
bined is about 16% of C2. The extremely large correction 
to ci can be simply explained. The cancellation which af- 
fects the single-particle contribution (see Eq. flBt ) does 
not apply to the two-body one. By using the DDH best 
values, all contributions due to PNC ECs give an asym- 
metry about 3.7xl0~ n , which is a 5% correction. The 
reason why large effects from individual meson exchanges 
lead to an overall small correction is due to the cancella- 
tion between pion and heavy-meson exchanges: the DDH 
best values have opposite signs for the pion and heavy- 
meson couplings. This conclusion however depends on 
the sign we assumed for the g p7r7 coupling. 
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Figure 4: The energy-dependences of PNC EC coefficients in 
the asymmetry parametrization: the top panel shows c 2 i in 
Eq. 11711 and the bottom panel shows c[ 2 2 3 4 6 of Eq. 1181 . 

IV. CONCLUSION 

The present work has been motivated by various as- 
pects of the PNC asymmetry A 7 in the deuteron photo- 
disintegration, especially in the few-MeV photon-energy 
range. A first work addressing this energy domain (l3l | 
showed that the process could provide information on the 
PNC ttNN coupling constant, h\, which allows one to 
check results from other processes involving this coupling. 
A later work rather schematic, concluded that this 
contribution could be largely suppressed. Between these 
two extreme limits, the question arises of what this con- 
tribution could be when a realistic description is made, 
including in particular the tensor and spin-orbit compo- 
nents of the NN interaction. At the first sight, a sizable 
PNC 7r-exchange contribution could arise if one assumes 
tensor-force effects of about 15% for each partial contri- 
bution and no cancellation. 

The complete calculation shows that the 7r-exchange 
contribution remains strongly suppressed after improving 
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upon the schematic model. Beyond making this obser- 
vation, a genuine explanation should therefore be found. 
When considering the asymmetry A 1 , an average is made 
over the spins of initial and final states. Terms in the 
interference effects of electric and magnetic transitions, 
whose spin dependence averages to a non-zero value, are 
expected to produce a sizable contribution. This discards 
the 7r-exchange contribution which involves a linear de- 
pendence on the spin operator S and tensor-force effects 
which involve the product of spin operators of order 1 
and 2. The argument applies to MECs too. A differ- 
ent conclusion would hold for an observable implying a 
spin polarization of the initial or final state. It thus ap- 
pears some similarity between the relative role of various 
contributions here and that one emphasized by Danilov 
for the inverse process at thermal energies: the circu- 
lar polarization of photons P 1 (equivalent to A 1 here) 
is mainly dependent on the PNC isoscalar and isotensor 
contributions while the asymmetry of the photon emis- 
sion with respect to the neutron polarization depends on 
the 7r-exchange contribution. 

As the 7r-exchange contribution to the asymmetry A 1 
turns out to have a minor role, we can concentrate on 
the vector-meson ones. At the low energies considered 
here, it is expected that these contributions depend on 
two combinations of parameters entering the description 
of the PNC (and PC) NN interaction. They are the 
zero-energy neutron-proton scattering amplitudes in the 
T = and T = 1 channels, At and A s . In terms of 
these quantities introduced by Danilov Q (see also later 
works by Missimer [2^], Desplanques and Missimer [2f| . 
Holstein j23) the discussion could be simpler. The asym- 
metry is found to vary between 

A~ ( = 0.70 m N X t — 0.17 m N X s at threshold 

and 

A 1 = -0.037 m N A t + 0.022 m N X s at uo 1 = 12 MeV, 

thus evidencing a change in sign which occurs around 
cj 7 = 5.5 MeV for both amplitudes. Depending on low- 
energy properties and, thus, on the best known properties 
of the strong interaction, the place where the cancellation 
of A 1 occurs sounds to be well established. It roughly 
agrees with what can be inferred from the analytic work 
by Khriplovich and Korkin ^4|. Not much sensitivity 
to PNC ECs is found. An experiment should therefore 
aim at a measurement at energies significantly different, 
either below or above. 

The goal for studying PNC effects is to get informa- 
tion on the hadronic physics entering the PNC NN in- 
teraction. This supposes that one can disentangle the 
different contributions to each process. We notice that 
the combination of parameters At and X s appearing in 
the expression of A 1 is orthogonal to that one deter- 
mining PNC effects in most other processes, especially 



in medium and heavy nuclei. The study of the present 
process is therefore quite useful. Another observation, 
which is not totally independent of the previous one, 
concerns the isotensor contribution. This one is espe- 
cially favored in the present process while it is gener- 
ally suppressed in processes involving a roughly equal 
number of protons and neutrons with either spin [2fl |. 
The present process is therefore among the best ones to 
get information on the isotensor pNN coupling constant. 
We however stress that this supposes the isoscalar parts 
could be constrained well by other processes. In a meson- 
exchange model of the PNC interaction, these ones are 
represented by the isoscalar pNN and uNN coupling 
constants. One could add that the relative sign of these 
two contributions is the same, in a large range of the 
photon energy (w 7 > 3 MeV), as in many other pro- 
cesses. It however differs at small photon energies where 
the asymmetry involves a combination of the various 
isoscalar and isotensor couplings that is little constrained 
by other processes. This explains that expectations of 
A 1 up to 10~ 7 near threshold could be suggested in re- 
cent works on the basis of a phenomenological analysis 
[lH. I28L l2fl|. Measuring this asymmetry could therefore 
be quite useful to determine a poorly known component 
of PNC NN interactions. On the theoretical side, the 
present work should be completed by the contribution of 
further parity-conserving exchange currents, but also by 
higher l/m^v-order corrections from the single-particle 
current and, consistently, from both PC and PNC ex- 
change currents Though they are not expected to 
change the main conclusions reached here, they could be 
required to obtain from experiment a more accurate in- 
formation on PNC NN forces. 
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Appendix: NON- VANISHING MATRIX 
ELEMENTS OF PNC ECS FOR DOMINANT 
TRANSITIONS 

In this section, we summarize the analytical expres- 

~ (s') ~ (2' ) 
sions of the non-zero F^^ and for the five domi- 

nant channels which lead to the numerical results in Sec- 
tion unci 
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1. F, 



bi 5 



As discussed in Section |n] an exchange charge at O(l) should contribute to this form factor. According to Ref. 
p3 |. the p-exchange does generate one: 



/r 
2V6 



Pm eso ™ c ( a; ; r i' r 2) = 2e 5pJVJV h° p --^= ] (ri x t 2 ) z ((Ti - <t 2 ) ■ V x [f p (r xl )f p (r x2 )) , (A.l) 



with f x {r) = exp(— m x r)/(Anr) and r xi = \x — r*|. The 1 S'o state is the only open exit channel and it gives 

(El? ) = 8 ^ ~ i^) ^ 5o l r /pMI 3 ^)^ ' ( A - 2 ) 
where (f\F(r)\i)d denotes the radial integral / drU*(f) F(r) Ud(i) and the subscript "<f refers to the deuteron state. 

As the four allowed exit channels 3 Pq, 3 Pi and 3 P 2 - 3 F 2 are spin- and isospin-triplet, the non-vanishing PNC ECs, 
which satisfy the spin and isospin selections rules, are 

j£rf>5»-i,r a ) = ^^hlf p (r)(r z -r z )(a 1 +a 2 )[(l + r^ 3 Hx-r 1 ))+(1^2), (A.3) 
^(xjn.ra) = ^^^/ w WW-r|)( < T 1 + < T 2 )((l + rn ( 5( 3 )( : r-r 1 ))+(1^2), (A.4) 

-Mv [(o-i x Vi) V 2 + <T?ff 2 x V 2 , / p (r s i) / p (r a2 )] ) + (1 <-» 2) , (A.5) 



3Zsoni c {x\r u r 2 ) = ^^/.i(T I xr 2 f(V 1 xV 2 )(/M,M) + (l-2). (A.6) 

V 2 ??i p v 7 

Note that one additional strong meson- nucleon coupling constant, g p7r7 , appears in Eq. l|A.6j) . This could be con- 
strained by the p — > it + 7 data. For the numerical calculation, we quote the number g p7rj = 0.585 as given in Ref. 

pill ]. The matrix element (Mlg 2 ') can be written as a sum of the contributions from each EC as 

(Mlf ] ) = — ( 5 plvjv ^ Ai + # WJV1V hi X 2 + g pNN ( h° p - ] A 3 ) + — 5 pJVJV g p7r7 h\ X 4 , (A. 7) 

"Ik y y 2v6 / / m p 

and for each exit channel, the quantities X\ 2 3 4 are 
1. 3 P 

*i = -|(( 3 Po|r/ P (r)| 3 5 1 ) d + i=( 3 Po|r/ p (r)| 3 ^ 1 ) d ) , (A.8) 
A 2 = ^( 3 P |r/.(r)| 3 5 1 ) (i + -i=( 3 Po|r/4r)| 3 J Di)^ , (A.9) 



Xa = 



■| ((1 + 2/V) ( 3 Pok /p(0| 3 5i), + -3 (1 - Mv) ( 3 P |r-/ P (r)| 3 A) d 

< 3 P |r/ P (r)| 3 S< +) ) d - ^ fPolr/pWI 3 ^"^ ] , (A.10) 

771p Tflp J 

{{ 3 P \f: p (r)\ 3 S 1 ) d -V2( 3 P \C(r)\ 3 D 1 ) d ) ; (A.ll) 



3(m 2 -ml) 



2. 3 Pi 
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X- 



3 



4 {( S Pi\rf p (r)\ 3 S 1 ) d -V2( 3 P 1 \rf p (r)\ s D 1 ) d ) , (A.12) 

^Wl^^d-x/lfPilr/^WI 3 ^)^ , (A.13) 

^((1 ~ Mv) ( 3 Pi|r / p (r)| 3 ^) d - V2(l - M ( 3 Pi |r / P (r)| 3 Di) d 

-A ( 3 Pl | r/p(r . ) |3 5 (+) )d + ^ ( 3 Pi|r/p(r)| 3^(-) ) \ (A _ 14) 



*4 = - , 



V3( 



((^il/JpMI 3 ^), + -L < 8 Pi|/;„(r)| 3 I>i> d ) ; 



3. 3 P 2 - 3 P 2 



X = 




x 2 = 






x 3 = 





P 2 |r / p (r)| 3 ^) d - ^ ( 3 P 2 |r/ p (r)| 3 ^) d - ^ ( 3 P 2 |r /.(r)! 3 ^}^ , (A.16) 
( 3 P2\rU(r)\ 3 S 1 ) d - ^(^rMrtfDiU - ^ ( 3 P 2 |r / w (r)| 3 J , (A.17) 



_^(l_ 4 ^)( 3 P 2 |r/ p (r)| 3 A) d 



1 3 , 



; < 3 P 2 |r /p (r)| 3 S«) d + M ( 3 P 2 |r/ p (r)| 3 p(-)) c 



?7ip 5 nip 



^(l + ^ v )^F 2 \rf p (r)\ 3 D 1 ) d + ^( 3 F 2 \rf p (r)\ 3 D^) d ) , (A.18) 



5 5 m p 



* 4 = 3(4^) [^'^^ iTf ( 3 P 2 |/; p MI 3 ^i) rf + ^ < 3 P 2 |/; p (r)| 3 A) d j , (A.19) 
where | 2S+1 L ( 7 +) ) = - ^) \ 2S+1 Lj), \^L^) = + | 2S+1 P/>, and /» ee f (/ T (r) - /»). 
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